Reconstructing f(R) theory according to holographic dark energy 
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In this paper a connection between the holographic dark energy model and the f(R) theory is 
established. We treat the f{R) theory as an effective description for the holographic dark energy 
and reconstruct the function f(R) with the parameter c > 1, c = 1 and c < 1, respectively. We show 
the distinctive behavior of each cases realized in f(R) theory, especially for the future evolution. 



I. INTRODUCTION 



energy is proposed by saturating the bound as 



Since 1998, the Type la Supernovae observations [lj 
have indicated that the expansion of the universe is cur- 
rently accelerating. This result has then been further 
confirmed by independent observations of Cosmic Mi- 
crowave Background (CMB)0 and Large Scale Struc- 
ture (LSS)[3(. One explanation for the cosmic accel- 
eration is ascribed to adding an exotic energy compo- 
nent with negative pressure, dubbed the dark energy, of 
which the origin and nature is still a mystery. Various 
dark energy models have been proposed in the litera- 
ture (see[J| for a detailed review). Among others, the 
simplest candidate for dark energy is the cosmological 
constant or the vacuum energy. Fitted quit well with ob- 
servational data though it is, this model suffers from the 
famous cosmological constant problem Q. This problem 
arises primarily due to the fact that the vacuum energy 
is considered within the framework of quantum field the- 
ory in Minkowski background. As we known, however, 
at cosmological scales where the effect of gravity has to 
be taken into account, the above description of the vac- 
uum energy would break down, and it is believed that 
the correct theoretical value of the vacuum energy will 
be predicted by a complete theory of quantum gravity. 
Although we are far from reaching such a fundamental 
theory, we do know some features of it. The holographic 
principle is an importance feature which can shed some 
light on the cosmological constant problem and the dark 
energy problem. According to this principle, considering 
gravity, the number of the degree of freedom of a local 
quantum field theory system is related to the area of its 
boundary, rather than the volume of the system as ex- 
pected when gravity is absent. Along this line, Cohen et 
al.Q suggested an entanglement relation between the IR 
and UV cut-offs due to the limitation set by the forma- 
tion of a black hole, which in effect sets an upper bound 
for the vacuum energy 
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where p\ is related to the UV cut-off, L is the IR cut-off 
and M p is the reduced Planck mass. The form of dark 
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where c is a numerical factor. It is easy to check that 
insert L = Hq 1 may give rise to a energy density com- 
patible with current observation in orders of magnitude. 
However, as Hsu[8j pointed out, this can not lead to a 
desired equation of state. Li0 proposed the holographic 
dark energy model, where L is chosen to be the future 
event horizon 
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This model has been tested to be well consistent with 
current observations [ldj . and it is a compelling candidate 
for solving the dark energy problem. 

Adding new component of dark energy to the whole 
energy budget is one way, there are also other promis- 
ing ways without resorting to new forms of energy. Since 
general relativity is only tested within solar system up 
to now, we may well consider the modification to the 
Einstein-Hilbert action in Einstein gravity at larger scales 
with higher order curvature invariant terms such as 
R 2 , R^R^u, R 1 ""* 13 Rpvap, or RD k R as well as nonmin- 
imally coupled scalar fields with terms like <p 2 R. Further- 
more, these terms naturally emerge as quantum correc- 
tions in the low energy effective action of quantum grav- 
ity or string theory [n], [T§]- Here we focus on the f(R) 
theories where the modification is a function of the Ricci 
scalar only. In[l3j|. the author first introduced an addi- 
tional R 2 to the Einstein-Hilbert action leading to an in- 
flationary solution for the early universe. As for applica- 
tions to the dark energy problem, it is shown that adding 
an term l/i? n [l4[ or more generalized Cx/R n + C2-R" 1 |l5l | 
can lead to late time acceleration originated from pure 
geometrical effect, equivalent to introducing an effective 
dark energy in the Einstein frame (see, for example, Sec. 
XVI. in [J] and the references therein for more works 
dedicated to solving the dark energy problem with f(R) 
theories). 

It can be shown flU] that through conformal trans- 
formations, extended theories with higher order terms 
and/or nonminimally coupled scalar fields correspond 
to Einstein gravity with some minimally coupled scalar 
fields (quintessence) suggesting, to some extent, an equiv- 
alence between dynamic dark energy models and f(R) 
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theories. An approach was proposed in[l7j to reconstruct 
the form of f{R) from a given expansion history H{z). 
From observational data such as SNe distance modulus 
vs redshift, we can obtain the luminosity distance Dl(z) 
through which H(z) is given by 
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In fact, however, large errors in current observational 
data prevent using this procedure to determine the exact 
form of f(R)- But the important thing is that we can 
use a given H(z) predicted from a dark energy model 
to reconstruct its equivalent f(R) theory. For example, 
the f(R) theories reconstructed from the H (z) given by 
quiessence and the Chaplygin gas respectively were ob- 
tained in[T3|. In this paper, we treat the holographic 
dark energy model as one inspired by the holographic 
principle, an important feature of a more fundamental 
theory of quantum gravity, and to reconstruct the corre- 
sponding f(R) theory as an equivalent description. Com- 
pared with previous works, where the holographic dark 
energy is reconstructed within scalar field models like 
ghost condensate [3, quintessence [l^, tachyon[2(| and 
hessence[2l|. here we perform the reconstruction in f(R) 
theory without resorting to any additional dark energy 
component, that is, the holographic dark energy is effec- 
tively described by the modification of gravity. 

In addition we note that there are in fact two strategies 
in f(R) theories: the metric formalism, where the action 
is varied with respect to the metric only; and the Pala- 
tini formalism [22|], where the metric and the connnection 
are treated as two independent variables with respect to 
which the action is varied. It is only in Einstein gravity 
f(R) = R that both approaches reach the same result. 
In general f(R) theories, the problem of which approach 
should be used is still an open question and the final so- 
lution may be determined by further observations and 
theoretical development. At present we assume the met- 
ric formalism in this paper. 



II. RECONSTRUCTION OF f(R) THEORY 

Now let's consider a homogeneous and isotropic uni- 
verse with flat spatial geometry consisting of matter and 
the vacuum energy given by @. The Friedmann equa- 
tion reads 
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where p m — p m o(l + z) 3 by the equation of energy con- 
servation, and a subscript denotes the value at present. 
By introducing fl A = 3 ^i H i and Sl m = 3 £? H i = 

^mo(l + z ) 3 j we obtain the Hubble parameter 



H(z) = H a 
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Clearly, once we determine the evolution of FIa (z), the 
whole expansion history H{z) is determined. Combining 
the definitions of holographic dark energy and the event 
horizon |(3J) and © we get 
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with the initial condition given by setting z = 1 in 
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Inserting ([6]) into the above equation and taking deriva- 
tive with respect to z on both sides(using 1 + z = 1/a), 
we obtain the differential equation of Oa 
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where the prime denotes derivative with respect to z. 
We can see that c is the only parameter affecting the 
dynamics of the holographic dark energy. In fact, we 
can use the equation of energy conservation of the dark 
energy to get the equation of state of the dark energy 
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It is easy to see that when — * 1 in the future, for c > 1 
the EoS will always be greater than —1 behaving like a 
quintessence; for c — 1 the universe will end up with a 
de Sitter phase; and for c < 1 the universe will end up 
with a phantom phase and the EoS crossing —1 occurs 
during the evolution exhibiting a quintom-like behavior. 
Therefore the parameter c plays a very important role in 
determining the evolutionary nature of the holographic 
dark energy. Many works [10j have been devoted to con- 
strain this parameter by observations such as SNe, CMB 
and galaxy clusters etc. Almost all the best fits indicate 
c < 1, although c > 1 is also compatible with the data 
within la. 

Once we fix c and fi m oj the evolution of Qa (z) can 
be determined by solving @ with the initial condition 
([5]). Then by © we can obtain H{z). Now we assume 
the holographic dark energy as an underlying theory of 
dark energy, and we want to find the corresponding f(R) 
theory as an effective description. We follow the method 
proposed in[T3|. In a FRW universe, the Ricci scalar can 
be expressed in terms of Hubble parameter: 
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In this paper we set k = 0. Note that here we assume 
the signature as {+,—,—,—}, the same as in[l7j. and 
therefore R is always negative (so it is with /). Once 
we choose { — , +, +, +}, the minus sign in front of (fTTT) 
would disappear. This is just a matter of convention, 
which means no physical difference. Let's start from the 
action 
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where C m is the matter Lagrangian. We use the units 
M p = c = H = 1. Variation with respect to the metric 
leads to the modified field equation [24j 
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T> rp^CU 



rp{m) 



(13) 



where G^y is the Einstein tensor, and an effective stress- 
energy tensor containing the higher order contributions 
is defined by 
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and the matter's contribution is in — T^ l> / f '(R) 

with the standard minimally coupled matter stress- 
energy tensor. With the FRW metric, we obtain the 
modified Friedmann equations 
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and the continuity equation 
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The three equations are not independent and we combine 
them to get one equation 
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where the prime denotes derivative with respect to R. 
Using the relation d/dt = —(1 + z)Hd/dz to replace the 
variable t by z , (fl8)) can be transformed into a third 
order differential equation of f(z) 

V . , (l9) 

where C n (z) consists of H (z) and its derivatives. Once we 
know the function H(z), the coefficients C n (z) , s can be 
calculated. However, in our case of the holographic dark 
energy model, H (z) cannot be derived analytically. By 
(I!]) and ([9|), H (z) and its derivatives can be expressed by 
JIa(-z), therefore, after painful derivation, the coefficients 
C n (z)'s can be ultimately expressed by the combinations 
of Qa(z). Although £!a itself can be solved by ([9]) alone, 
it has to be considered as a part of the whole set of dif- 
ferential equations in order to solve f(z). So we consider 
© and (TIT?|) as a differential equation set, of which one 



initial condition is ([8]) . According to [l7| , the other initial 
conditions are 
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f(z = 0) = f(R ) = 6ff 2 (l - n m0 ) + R . 
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Given Jl m o and c, with these four initial conditions 
©and ((101)- d22]), the differential equation set ® and 
(112]) can be solved numerically. The reconstructed func- 
tion ,f(R) is presented in FigfTl where we set f2 m o = 0.29 
and c = 0.6, 0.8, 1.0, 1.2 respectively. For the sake of 
comparison, we also show the same result on a If — IR 
plane in Fig(2l where If = ln(—f) and IR = ln(—R) as 
used in[T3|. Note that we set the values of c within the 
range 0.6 < z < 1.2, which is consistent with fitting re- 
sults according to the works in[10j. Compared with the 
reconstructed function f(R) for quicssencc and that for 
Chaplygin gas (Fig.2 and Fig.4 in[I3), we can see that 
the three figures are similar. This is because the values 
of the parameters are set to be around their best fits. 
This is in effect approximately equivalent to reconstruct 
f(R) with H(z) directly from observational data. We 
expect future observations with more accurate data will 
discriminate between these models. 
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FIG. 1: Reconstructed f(R) with < z < 10 and c = 
(dash-dotted), c = 0.8 (dahsed), c = 1.0 (solid) and c 
(dotted). 
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III. DISCUSSION AND CONCLUSION 

The reconstructed f(R) theory is naturally consistent 
with the solar system experiment by the reconstruction 
method. In fact, this requirement is just the physical 
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FIG. 2: Reconstructed f(R) in If — IR plane with < z < 10. 
If = ln{-f) and IR = ln{-R). 



motivation for the initial conditions (j2"0|) and (j2"Tj) . 
one thing, if we rewrite (|15[) explicitly with 87rG 



if 2 



8nG 



/'(£) 



For 



(23) 



we can see that f'(R) is effectively modified the Newton 
gravitational constant G as G/f'(R), that is , a vari- 
able gravitational coupling. In order to be compatible 
with solar system experiments, at z = 0, we must require 
G/f(R ) = G or f(R a ) = 1. By 



f'(Ro) = i 



dz 



= 1. 



(24) 
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this leads to (|2U)l . For the other, it is shown (2o| that the 
consistency with solar system test also requires f"(Ro) = 
0, which directly gives rise to (f2Tj) . As for the remote 
past, there is no reason for us to impose this requirement 
since the experiments are done today and the validity of 
the result holds only at z ~ 0. 

FigQ] shows that for small \R\ (small z also), the func- 
tions f(R) are indistinguishable for different parameter c. 
As we mentioned before, c is a crucial parameter charac- 
terizing the nature of the holographic dark energy. Differ- 
ences between the corresponding f(R) functions become 
significant as \R\ (or z) increases. To further illustrate 
that the reconstructed theory does reflect the distinc- 
tive effect of c, we consider the future evolution scenario. 
FigEl shows the future evolution of R. As is expected, 
for c < 1, the curves indicate the typical phantom be- 
havior: \R\ — > oo. This is because the dark energy with 
EoS< —I dominates over matter and the phantom en- 
ergy density increases with time, tears apart structures 
and a Big Rip is unavoidable. For c = 1, \R\ varies little 
and the dark energy becomes more and more like a cos- 
mological constant. For c > 1, R vanishes in the future. 
In FiglU we can see that the difference is more distinc- 
tively reflected by the function f(R) reconstructed based 



upon the future evolution of the holographic dark energy 
model. For c = 1.2, as R approaches zero, / increases 
from negative to positive, which may indicate a inverse 
power law dependence of f on R. This is consistent with 
the models proposed in [14[ and [HI]. For c = 1, the 
straight line manifests a linear dependence on R up to 
a constant, which is consistent with the de Sitter phase 
where f = R + 2A. For c = 0.8, the curve first meets 
a turnaround point, at which the decreasing \R\ begins 
to increase due to the domination of the phantom-like 
dark energy. It can be checked that in this case, the 
turnaround redshift is in the near future for c = 0.8 while 
it is in the near past for c = 0.6, namely, the domination 
of phantom-like dark energy begins earlier for smaller c. 
As the universe evolves, \R\ keeps growing, and / de- 
crease first and then increases to become positive. Both 
/ and R become divergent in the final Big Rip. Further 
analysis shows that for c = 0.6, / — > — oo; for c = 0.8, 
/ — > +oo, which clearly reveals that the form of the func- 
tion f(R) may be significantly different for different c. 
We note that the existence of the turnaround point is a 
universal feature for all the phantom-dark energy models 
realized in f(R) theories, due to the competition between 
dark energy and matter. 
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FIG. 3: The future evolution of R 



In conclusion, we have reconstructed the function f(R) 
in the extended theory of gravity according to the holo- 
graphic dark energy. The basic reconstruction procedure 
can be simply summarized as: first, express R as a func- 
tion of z\ then / can also be considered as a function 
of z by f(z) = f[R(z)); thirdly obtain a third order dif- 
ferential equation for f(z) and solve it with some initial 
conditions; and finally reconstruct f(R) from R{z) and 
f(z). Note that in this procedure, the Hubble parameter 
H(z) and its derivatives with respect to z enter into R 
and the coefficients of the differential equation for f(z). 
Since H(z) and its derivatives can be expressed by 
what we are dealing with can be treated essentially as a 
differential equation set with the unknown functions f{z) 
and £l\(z) to be solved. With some initial conditions im- 
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FIG. 4: Reconstructed f(R)- The curves are plotted with z 
from around 2 down to —1. The arrow denotes the direction 
of 2 increasing. The star denotes current value at z = 0. 
The turnaround point is the point where the total energy 
density starts to increase after phantom energy dominates 
over matter. 



posed by physical consideration, we solved the differential 
equation set and found the function f(R) numerically. 

In addition, it should be emphasized that the holo- 
graphic dark energy is the result obtained within the 
framework of general relativity, rather than any other 
extended theory such as f(R) theory. What we have 
done is to reconstruct the f{R) theory which effectively 



describes the holographic dark energy in Einstein gravity. 
Whether the holographic vacuum energy can be general- 
ized to f(R) theories and what it looks like are questions 
worth further investigation. 
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APPENDIX 

Here we list the some expressions essential for practical 
calculation. The coefficients in (1191) are 
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R(z) and its derivatives are 
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Higher order derivatives of H(z) are too complicated to 
be listed here. In practical calculation, we use com- 
puter program for derivation. In addition, the code for 
Mathematica 5.0 we used is available on request. 
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